Ground-state properties of a one-dimensional system of dipoles 
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A one-dimensional (ID) Bose system with dipole-dipole repulsion is studied at zero temperature 
by means of a Quantum Monte Carlo method. It is shown that in the limit of small linear density the 
bosonic system of dipole moments acquires many properties of a system of non-interacting fermions. 
At larger linear densities a Variational Monte Carlo calculation suggests a crossover from a liquid- 
like to a solid-like state. The system is superfluid on the liquid-like side of the crossover and is 
normal in the deep on the solid-like side. Energy and structural functions are presented for a wide 
range of densities. Possible realizations of the model are ID Bose atom systems with permanent 
dipoles or dipoles induced by static field or resonance radiation, or indirect excitons in coupled 
quantum wires, etc. We propose parameters of a possible experiment and discuss manifestations of 
the zero-temperature quantum crossover. 
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Up until now Bose-Einstein condensation has been re- 
alized in many different atom and molecule species with 
short-range interactions. At low temperatures such an in- 
teraction can be described by a s-wave scattering length 
and is commonly approximated by a contact pseudopo- 
tential. In contrast, some recent work has focused on the 
realization of dipole condensates 0, El S 01 • In these 
systems, the dipole-dipole interaction extends to much 
larger distances and significant differences in the proper- 
ties (such as the phase diagram and correlation functions) 
are expected. Another appealing aspect of a system of 
dipole moments is the relative ease of tuning the effec- 
tive strength of interactions which makes the system 
highly controllable. Dipole particles are also considered 
to be a promising candidate for the implementation of 
quantum-computing schemes 0, IE E} • 

On the theoretical side dipole condensates have been 
mainly studied on a semiclassical (Gross-Pitaevskii) [ljj 
or Bogoliubov level. A model Bose-Hubbard Hamil- 
tonian has been used to describe a dipole gas in optical 
lattices and a rich phase diagram was found [5|, [l2| ■ So 
far there have been no full quantum microscopic compu- 
tations of the properties of a homogeneous system. 

Recently the Monte Carlo method was used to study 
helium and molecular hydrogen in nanotubes|l.lj. Such 
a geometry, which is effectively one dimensional, leads 
to completely different properties compared to a three- 
dimensional sample. 

We consider N repulsive dipole moments of mass M 
located on a line. The Hamiltonian of such a system is 
given by 
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We keep in mind two different possible realizations: 
1) Cold bosonic atoms, with induced or static dipole 
momenta, in a transverse trap so tight that excitations 
of the levels of the transverse confinement are not possi- 
ble and the system is dynamically one-dimensional. The 



dipoles themselves can be either induced or permanent. 
In the case of dipoles induced by an electric field E the 
coupling has the form Cdd = E 2 a 2 , where a is the static 
polarizability. For permanent magnetic dipoles aligned 
by an external magnetic field one has Cdd = rn 2 , where m 
is the magnetic dipole moment. The interaction between 
the atoms, apart from the dipole forces, contains a short- 
range scattering part which is conveniently described by 
the s-wave scattering length a. Usually this contribution 
is large compared to the dipole force, but recent progress 
in applying Feshbach resonance techniques to tune o, and 
even to make it zero, opens exciting prospects of obtain- 
ing a system with purely dipole-dipole interactions. The 
strength of the effective coupling Cdd can be tuned by 
changing the electric field in the case of induced dipoles, 
and the special technique of fast rotation of the electric 
or magnetic field can be applied to permanent dipolesQ. 

2) Spatially indirect excitons in two coupled quantum 
wires. A quantum wire is a semiconductor nanostruc- 
ture where an electron or a hole is allowed to move only 
in one direction and excitations of the transverse quan- 
tization levels are negligible. In two parallel quantum 
wires, one containing only holes, and the other only elec- 
trons, holes and electrons couple forming indirect exci- 
tons. If such a system is dilute enough, it constitutes a 
ID set of dipoles. In this case Cdd = e 2 D 2 /e, where e 
is an electron's charge, e is the dielectric constant of the 
semiconductor, and D is the distance between the cen- 
ters of the quantum wires. This system is ID counter- 
part of 2D indirect exciton system in coupled quantum 
wells, which was extensively studied both theoretically 
[3 El 03 E3 I3 and experimentally 0. The prop- 
erties of ID and 2D systems differ essentially (see be- 
low Tonks- Girardeau regime, etc.), thus an experimental 
study of the ID system is of a great importance. 

The Hamiltonian Q can be written in dimensionless 
form by expressing all distances z = z/rQ in units of 
r = MCddj (27rfi, 2 ) and energies in units of £q — Ti 2 /Mr 2 . 
All properties of such a system are defined by the dimen- 
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sionless parameter nr$ with n — N/L being the linear 
density and L the size of the system. 

Our aim is to determine the ground state energy in 
wide range of densities and to check the evidence of a 
possible quantum crossover. In order to define the struc- 
tural properties we calculate the pair distribution func- 
tion (PDF). In terms of the many-body ground state 
wave function of the system \l/o PDF is written as 
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The static structure factor is directly related to the 
pair distribution function 



S(k)=l + n / e tkz {g 2 (z) - l]dz 



(3) 



The technique of Bragg scattering provides access to 
the static structure factor in experiments on bosonic 
atoms. For the system of indirect excitons measurements 
of spatial structure of photoluminescence can give infor- 
mation about possible crystallization, etc. 

We apply Diffusion Monte Carlo (DMC) method, 
which is one of the most efficient theoretical tools for 
investigating zero temperature properties of quantum 
systems 0|- We choose the many-body guiding trial 
wave function in the Bijl-Jastrow form consisting of 
one- and two- body terms: 



N 



*r(zi,. 



ZN) 



Zk 



(4) 



j<k 



To describe the liquid-like side of the crossover it is suf- 
ficient to have only the term describing pair-correlations. 
We choose it in the form /2(|z|) = exp{ — L4/(n|z|)] B }, 
where A and B are variational parameters which we op- 
timize by minimizing variational energy by carrying out 
Variational Monte Carlo (VMC) calculation. As will 
be discussed below, in the low density limit the wave 
function asymptotically approaches the Tonks-Girardeau 
gas wavefunction / 2 TG (|z|) = \sm{-Kz/L)\ [13, EH- We 
checked that in this regime there is no large difference 
between using /J G (|z|) and /zd^l) even on a variational 
level. Thus our choice of a two-body term is well suited to 
describing a liquid even in the strong mean-field regime 
(see IHl). 

We account for the spatial quasi-crystalline order by 
considering Gaussians with width C for each particle near 
to a corresponding lattice site = exp{ — [n(zi — 

zf)] 2 / (2C 2 )} . The sites zf are equally spaced with the 
distance n^ 1 and the variational parameter C is defined 
by VMC optimization. For the simulation at high density 
we keep the same type of two-body term f 2 as on the 
liquid-like side of the crossover, although the values of 
optimal parameters may differ. 

In Fig. we show the comparison of the VMC en- 
ergy of the liquid-like and solid-like wavefunctions in the 
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FIG. 1: Difference between energy calculated using liquid-like 
and solid-like wavefunctions (E so nd — Ei iqu id) / E as a function 
of the dimensionless density. 



nro = 0.01 — 0.1 density range. In this region the op- 
timal parameters are A — 1.6, B — 0.4, C = 1.16. We 



discover that for densities smaller than 



0.085rn 



the liquid-like wavefunction description energetically fa- 
vorable, although at larger densities the ground state of 
the system is better described by a solid-like ansatz. On 
the contrary to 3D and 2D liquid-solid phase transitions, 
the energetic difference in ID system is very small. This 
suggests that the transition in a one-dimensional system 
at zero temperature is of a crossover type. 

At a very small density one expects that: (1) in the 
process of two-body collisions particles always get re- 
flected back due to the repulsive interaction which is very 
intense at short distances, (2) particles stay far apart 
most of the time, so that the potential energy of the in- 
teraction is small compared to the kinetic energy. From 
above it follows that the system is equivalent to a gas 
of impenetrable bosons (Tonks-Girardeau gas). It was 
shown in |20j that the wave function of such a Bose sys- 
tem can be mapped onto the wave function of a system 
of non-interacting spinless fermions. The bosonic system 
acquires many fermion-like properties (fermionization) : 
the energy is the Fermi energy E TG /N = £oir 2 (nro) 2 /6, 
and the pair distribution function exhibits Friedel-likc 
oscillations gJ G (z) = 1 — sin 2 (nnz) / (wnz) 2 . 

We present the dependence of the energy per particle 
on the density in Fig.Q] At small density nro <c 1 the en- 
ergy is the same as that of the Tonks-Girardeau gas E TG , 
which signals fermionization of the system in this regime. 
At large density the particles are localized at lattice sites 
with the potential energy E strAnt '/N = £ ((3){nr ) 3 be- 
ing dominant. In this regime the density dependence of 
the energy is very strong. Indeed it is cubic, in contrast to 
the linear dependence on the mean field for the pseudopo- 
tential interaction and the quadratic dependence for an 
exactly solvable model with a 1/z 2 interaction [23|. This 
strong dependence comes from the diverging short-range 
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FIG. 2: Energy per particle as a function of the dimension- 
less density (solid line), energy E TG of the Tonks-Girardeau 
limit (dashed line), potential energy in the strongly interact- 
ing limit (dash-dotted line). Everything is measured in units 
of h 2 /mr 2 . 
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FIG. 3: Pair distribution function obtained from a DMC 
calculation for densities nro = 10 -3 ,0.1, 1, 10 (larger densities 
have higher peaks). 



behavior of the dipole-dipole interaction. Although in 
3D the dipole-dipole interaction is long-range, this is no 
long er true in the one-dimensional case (see, for example, 
|3) and no special techniques like Ewald summation are 
required. 

We studied the dependence of the energy on the size 
of the system. The energy has two contributions: one 
coming from summation over pairs separated by a dis- 
tance smaller than L/2 (this contribution is a direct out- 
put of the Monte Carlo calculation) and a tail energy, 
which is estimated by approximating the density at L/2 
by the asymptotic bulk value. We find that the energy 
per particle as a function of the system size quickly satu- 
rates to its thermodynamic value and the results obtained 
for N = 50, 100, 200 particles agree within the statistical 
accuracy present in our calculation. We also find that 
the "tail energy" contribution is smaller than 0.2% of 
the total energy. All reported results are obtained using 
N = 100 particles. 

The pair distribution function (Eq. (0)) is presented in 
Fig. |31 for a range of densities covering all the regions of 
interest. In a dilute system we find that amplitude of the 
oscillation decays rapidly which is characteristic for a liq- 
uid. In particular, at the density nro = 10~ 3 it is almost 
impossible to distinguish the pair distribution function 
from that of spinless fermions, which is in agreement with 
the arguments given above. By increasing the density the 
oscillation becomes more pronounced. Around the criti- 
cal density n c r w 0.08 we enter the solid-like side of the 
crossover. Further, at larger densities nro = 1, 10 we see 
manifestations of the localization of particles near lattice 
sites. 

By performing the Fourier transformation Q we ob- 
tain the static structure factor from the pair-distribution 
function. Knowledge of S(k) is of high importance as 
it can be accessed experimentally by using Bragg spec- 
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FIG. 4: Static structure factor obtained from a DMC calcu- 
lation for densities nro = 10 -3 , 0.1, 1. A higher first peak 
corresponds to a higher density. 



trometry. At small densities the low-momentum region 
has a linear form which is related to the speed of sound c 
by the Feynman relation S(k) — h\k\/(2Mc). In par- 
ticular, in the Tonks-Girardeau limit the static struc- 
ture factor takes an extremely simple form: the linear 
growth matches the asymptotic constant at the wave 
vector |fc| = 27m. Increasing the density leads to the 
formation of a peak structure. Presence of the peak in 
the S(k) is a consequence of the dipole-dipole interac- 
tion, as the peak is absent in a system with only s-wave 
scattering |22| . In the regime of large density we see sev- 
eral peaks at integer multiples of 27m. 

In order to test superfluidity of the system we use 
the winding- number technique|25j. The superfluid frac- 
tion is given by the ratio of the imaginary time diffu- 
sion constants of the center of mass and of a free par- 
ticle [2(| . In contrast to the calculation of the energy, 
here the symmetry of the trial wavefunction is crucial. 
We restore symmetry in the trial wavefunction in the 
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solid-like region by performing a summation over all sites 
= Ei exp{-[n(^ - z|)] 2 /(2C 2 )}. We see negligi- 
ble differences in the energy (similarly to that found in 
4 He simulations which justifies our previous choice 
of ^>t- We find that the system is superfluid on the 
liquid-like side of the crossover and is normal deep on 
the solid-like side. 

We argue that the critical densities of the quantum 
crossover can be readily reached in experiments with in- 
teracting indirect excitons in two quantum wires. Let 
us take the GaAs parameters (e = 12.5, electronic mass 
m e = 0.07mo, hole mass mh = 0.15too, with mo being 
the free electron's mass) for reference. Then the mass of 
an exciton is 0.22mo. Using quantum wires with a sep- 
aration between their axes, D, equal to 5nm, we obtain 
ro ~ 10~ 6 cm and £q w 3meV\ Then the dimensionlcss 
density nro = 0.08 — 0.09 corresponds to an achievable 
linear excitonic density n rs 10 5 cm _1 . The liquid-like 
and solid-like regions and the crossover between them 
should be realized by changing the density of excitons in 
the quantum wires. 

Interactions between chromium atoms can be effi- 
ciently tuned as proposed in [||. Chromium has the ad- 
vantage of having a large permanent magnetic moment 
when compared to the other alkali atoms. Research into 
achieving Bose-condensation in chromium is now a hot 
topic P, 0, El [1- I n permanent dipole chromium atoms 
the ratio between the strength Cdd of the dipole-dipole 
interaction and s-wave coupling constant is 0.27 for 52 Cr. 
Manipulation of induced electric dipoles is more difficult, 
although in such a system a ratio of the order of 10 2 can 
be reached leading to the realization of an almost pure 
dipole system. 

In conclusion we have investigated the ground state 
properties of a system of dipole moments by means of 
a quantum Monte Carlo method. We have found the 
presence of a quantum crossover: at small linear densi- 
ties nro system stays liquid-like and superfluid, although 
in a more compressed system the solid-like description 
is energetically favorable. We have calculated the exper- 
imentally accessible pair distribution function and the 
static structure factor for a wide range of densities. In 
the dilute limit this Bosonic system becomes similar to 
a system of spinless fermions (fermionization) and the 
properties of the system are those of a Tonks-Girardeau 
gas. On the liquid-like side of the crossover the system is 
superfluid, although it is normal deep on the solid side. 

Finally, we have pointed out that the critical density of 
the quantum crossover can be achieved in current exper- 
iments with interacting excitons in wires and have pro- 
posed the parameters of a possible experimental setup. 
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